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We consider the evolution of a family of two-dimensional (2D) dispersive turbulence mod- 
els. The members of this family involve the nonlinear advection of a dynamically active scalar 
field, and as per convention, the locality of the streamfunction-scalar relation is denoted by 
a, with smaller a implying increased locality (a = 1 gives traditional 2D dynamics). The 
dispersive nature arises via a linear term whose strength, after non-dimensionalization, is 
characterized by a parameter e. Setting < e < 1, we investigate the interplay of advection 
and dispersion for differing degrees of locality. Specifically, we study the forward (inverse) 
transfer of enstrophy (energy) under large-scale (small-scale) random forcing along with the 
geometry of the scalar field. Straightforward arguments suggest that for small a the scalar 
field should consist of progressively larger isotropic eddies, while for large a the scalar field 
is expected to have a filamentary structure resulting from a stretch and fold mechanism; 
much like that of a small-scale passive field when advected by a large-scale smooth flow. 
Confirming this, we proceed to forced/dissipative dispersive numerical experiments under 
weakly non-local to local conditions (i.e. a < 1). In all cases we see the establishment of 
well-defined spectral scaling regimes. For e ~ 1, there is quantitative agreement between 
non-dispersive estimates and observed slopes in the inverse energy transfer regime. On the 
other hand, forward enstrophy transfer regime always yields slopes that are significantly 
steeper than the corresponding non-dispersive estimate. At present resolution, additional 
simulations show the scaling in the inverse regime to be sensitive to the strength of the 
dispersive term : specifically, as e decreases, quite expectedly the inertial-range shortens but 
we also observe that the slope of the power-law decreases. On the other hand, for the same 
range of e values, the forward regime scaling is observed to be fairly universal. 

PACS numbers: 47.52. +j 

I. INTRODUCTION 



It is quite common in geophysical fluid dynamics to encounter problems that involve the 
presence of both advection and dispersion (see for example, Chapter 5 in the text by Majda 
PQ and Chapter 5 in the text by Chemin et al. [2J). In a two-dimensional (2D) context, simple 
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model equations that possess advection and dispersion include the familiar barotropic beta plane 
equations [3] and the dispersive surface quasigeostrophic (SQG) equations [4J. An interesting 
aspect of such equations, highlighted by Rhines [3J in the barotropic beta plane case, is the 
co-existence of turbulent and wave-like motions [5], [6j. In fact, the dispersive SQG equations have 
been proposed as an alternate (more local) platform for investigating wave-turbulence interactions 

in- 



In the present work, we look at an extended family (that includes the aforementioned examples 
as members) of dispersive active scalars. This is a simple dispersive generalization of the family 
of 2D turbulence models introduced in Pierrehumbert, Held & Swanson [7J. Specifically, in a 2D 
periodic setting we consider 
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where 9{k x ,k y ,t) = -(k z ) a ^(k x ,k y ,t) (1) 

where D/Dt denotes the 2D material derivative and e is a non-dimensional parameter. Note 
that in real space 9 and ip above are related via a suitable (pseudo-) differential operator, i.e. 
9 = — (— A) a ip where A is the 2D Laplacian, and for our purposes a E For the beta 

plane equations we have a = 1, is the vorticity field, and e corresponds to the Rhines number 
Rh = U/{j3L 2 ). In the dispersive SQG case a = 1/2, 9 represents the buoyancy (or potential 
temperature [4J), and e = U/(AL). Physically, of course we are dealing with very different 
scenarios wherein (3 is the ambient planetary gradient of the vorticity while A is the background 
surface buoyancy gradient. Apart from a = 1, 1/2, there exist other members of this family that 
have physical interpretations [8J. But, from a broader perspective, the entire family is well defined, 
and we expect that varying a would provide greater insight into the interplay between advection 
and dispersion. 



We begin examining the effect of a on the dynamics by suppressing the dispersive term. After 
developing a feel for the dependence of the energy/enstrophy transfer and the geometry of the 
scalar field on a, we introduce dispersion and note its modulating effect on the nonlinear term. 
In the limit e — > 0, the so-called limiting dynamics are controlled by resonant interactions. In 
fact, the resonant skeleton corresponding to ([!]) is seen to be of the same form as that of the 
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rotating three-dimensional (3D) Euler equations. Following the original argument by Rhines [3], 
for < e < 1, if a > 0.5, the simultaneous transfer of energy to large scales and small frequencies 
leads to an anisotropic streamfunction (resulting in predominantly zonal flow). In contrast, given 
the nature of the dispersion relation, it turns out that for a < 0.5 the constraints of energy transfer 
to large scales and small frequencies do not require the spontaneous development of anisotropy in 
the streamfunction. Confirming this via decaying runs from spatially un-correlated initial data, we 
proceed to a suite of forced/dissipative numerical simulations. In all cases, under weakly-local to 
local conditions (i.e. a < 1) we note the formation of clear power-law scaling. The various slopes 
are extracted and compared with non-dispersive estimates. The paper ends with a collection of 
results and a brief discussion of potentially interesting avenues for future work. 



II. SOME OBSERVATIONS ON THE EFFECT OF a 



Before considering the effect of dispersion, we examine the influence of a in greater detail. 
To get a feel for the locality of the interactions, we shut off the linear dispersive term in 
and consider a simple numerical exercise involving the evolution of a smooth 9 ring. As is 
seen in Fig. 0, the deformation of the ring is more physically local for smaller a. Given that 
i/j(k) = —9(k)/k 2a , as a increases, one sees that only the small k features of 9 and i/j remain 
dynamically active. As a result, for larger a:, smaller scale features of the scalar field are for all 
purposes driven in a passive manner. In fact, in the limit a — > oc, only the wavenumber one 
{k — 1) component of 9 and i(j is coupled, and we end up with the problem of passive advection 
via a large-scale smooth flow. In contrast, the case of small a is markedly different. Specifically, 
(u,v) = (ik y 9/k 2a ,-ik x 9/k 2a ) results in a transition at a = 1/2 (when the scalar and velocity 
fields have similar scales), where for a < 1/2, the velocity fields are in fact of a finer scale than the 
advected scalar field. Note that ip remains a smoothened form of 9. We expect this transition to 
have interesting consequences as there are well-known differences in the behavior of a passive field 
when advected by large-scale, comparable-scale and small-scale flows respectively [9].[T0].fTT]. 

In a 2D periodic domain, much like its non-dispersive counterpart, ([!]) conserves energy 
E — J D ^9 and enstrophy or v 9- variance" £ — J- D 9 2 . Continuing with the non-dispersive case, 
following arguments for the 2D Euler equations [12J, it is expected that the energy primarily 
flows to large scales while the ^-variance is transferred to small scales. Assuming the presence of 
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FIG. 1: Non-dispersive evolution of rings. From left to right, a = 0.5, 1 and 1.5 respectively. Quite 
clearly, for smaller a the deformation is more local in character. 

an equilibrium cascade, detailed spectral scaling laws in the appropriate inertial-range have also 
been derived [7], jS], Q2] (see also Tran [H] and Tran & Bowman [15] for bounds on the scaling 
exponents for certain special dissipative operators in the SQG case). If we consider the Fjortoft 
estimate ([IB]. Chapter 4 in Salmon [17] ) i.e. the transfer of energy out of scale k\ into scales feg? ^2 
(s.t. fco = k\/2a and &2 = 2afci, where a > 1), we have Eq = E\ x {a 2a /[l + a 2a ]}; i.e. the more 
nonlocal the situation (a > 0.5), the larger is the fraction of the energy (enstrophy) transferred 
to large (small) scales. Note that increasing a shows this unequal distribution of energy and 
enstrophy is further exaggerated when the exchange involves scales that are progressively further 
apart, i.e. the energy /enstrophy partition is more biased in spectrally nonlocal transfers. 

To develop a feel for the dependence of the geometry of an emergent scalar field on a, we 
continue with non-dispersive simulations, though now from spatially un-correlated initial data 
chosen from a Gaussian distribution with unit variance. Given the presence of an inverse transfer 
of energy, we expect coherent structures to emerge from this un-correlated initial condition. In 
fact, given the similar scales of the velocity and scalar fields for smaller a we do not expect the 
scalar field to undergo much stretching and folding, while for large a we expect repeated events 
of this sort leading to a filamentary scalar field — much like the fate of a small-scale passive blob 
when advected by a large-scale smooth flow — due to the implicit large-scale strain (via the large 
scale-separation between 9 and ^). These expectations are confirmed in Fig. ^ which shows the 
initial condition and the emergent scalar field for a = 0.5 and 2 respectively. 
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FIG. 2: The first panel is the spatially un-correlated initial condition (smoothened via a diffusive stencil). 
The second and third panels show the emergent scalar field for a = 0.5 and 2 respectively. Quite clearly, 
for a = 0.5 we have a field composed of coherent eddies while for a = 2 we obtain a filamentary geometry 
reminiscent of a passive field when subjected to large-scale advection. 

III. THE INFLUENCE OF THE DISPERSIVE TERM 

When the dispersive term is included, the linearized form of ([!]) supports waves with the dis- 
persion relation 

k 

u ( k ) ; where k i k x,ky) (2) 

It is interesting to note that, for a = 0.5, ^ gives < \w(k)\ < 1/e. For a < 0.5, \u\ — > oc for 
k y = 0, k x — > oo, whereas for a > 0.5, \uj\ — > oc for k y = 0, k x — > 0. The dispersion relation for 
different a (on either side of a = 0.5) is plotted Fig. Q. Quite clearly the frequencies have a 
very different dependence on wavenumber when a <, = or > 0.5. In fact, an important feature of 
the beta plane equations (which is true for all a > 0.5) that |cj(fe)| increases for large scales is no 
longer true for a < 0.5; in fact, for a < 0.5 the smaller scale features have larger frequencies. 

In the following, we occasionally refer to modes with zero frequency as slow modes and those 
for which oo ^ as fast modes. In a periodic setting, substituting a Fourier representation in 
yields 

9 $k ii r _i r Px Qx . k x 

~0f =2^ C kp q Mp ex p{ — ["^ " ^ + p^]*} (3) 
where J2t represents a sum over q such that k = p '+ q and Ck pq is the interaction coefficient 
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FIG. 3: The dispersion relation ^ with e = 1 for, from left to right, a = 0.25, 0.5 and 1 respectively. Note 
that for a = 0.5 the frequencies are bounded while on either side we obtain — > oo in particular limits. 



given by 



Ckpq — (PxQy — PyQx)[ 



p 2a _ q 2a 



(4) 



A. Zonal flows 



Considering an initial value problem, the original deduction of anisotropic fields by Rhines [3J, 
in the context of the beta plane equations was based on the dual constraint of an upscale transfer 
of energy along with the tendency of resonant triads to move energy into small frequencies (see 
Hasselmann [18J for a general demonstration irrespective of the details of the nonlinear coupling co- 
efficient). Indeed, the two pieces of the Rhines argument are : (i) energy moving to large scales as a 
result energy/enstrophy conservation, and (ii) the importance of resonant triads in energy transfer 
when dispersion modulates the advective nonlinearity. In the context of the general dispersion 
relation pj), for isotropic structures, k x ~ k y ~ k we have 



Hence for a > 0.5, moving to large scales, i.e. for decreasing k we encounter larger frequencies. 
Therefore, to satisfy the dual constraints, Rhines [3] suggested that the system would sponta- 
neously generate anisotropic structures; further examining ^ shows that these constraints are 
satisfied for k y ^ Q,k x /k y 1. Also, when considering energy transfer into large scales, i.e. 
k < p,q, interactions that fall in the aforementioned anisotropic category are in fact near-resonant 
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[15] . In essence we have an anisotropic streamfunction that results in predominantly zonal 
flows. Though note that for a < 0.5, decreasing k implies smaller frequencies, therefore it is 
possible to maintain isotropy while simultaneously transferring energy to large scales and small 
frequencies. Hence, for a < 0.5, the dual constraints do not nessecitate the formation of a 
dominant of zonal flow. Note that this does not imply zonal flows cannot form for a < 0.5, in 
fact, substituting an expansion of the form = ^° + e^ 1 ..., the 0(1/ e) balance in yields 
dtp°/dx = =>* u° = u°(y,t),v° = 0. Of course, this expansion doesn't imply any control over the 
higher order terms, but irrespective of a, at order zero, it indicates the possibility of zonal flow 
formation. 

A different, though complementary, approach to the generation of zonal flows is to view as 
a mixing problem, i.e. to consider the mixing of 9 in the presence of a background gradient |20j . 
In the context of the beta plane equations, as shown by Rhines & Young [20], steady flows with 
closed streamlines (or gyres) result in the homogenization of potential vorticity between these 
streamlines. With a linear background, this naturally leads to a "saw-tooth" vorticity profile with 
high gradients concentrated on the streamlines [2T] . Further in a statistically steady state, this 
reasoning posits the existence of successive "saw-tooth" structures in the vorticity profile (or in 
other words a potential vorticity "staircase") if the size of the eddies is smaller than the size of 
the domain [21] (see [22], [23] for forced- dissipative simulations and [24J for experimental results 
in this regard). From this viewpoint localized zonal flows (jets) result from an inversion of the 
aforementioned potential vorticity "staircase". For the general family in 0, 9 + y/e is mixed by 
the flow and a similar homogenization within streamlines is expected to follow. Further, as the 
scale of the velocity field becomes smaller than that of the advected scalar when a crosses 0.5 
from above, we expect the homogenization for small a to progress in a more diffusive manner, 
i.e. small-scale inhomogeneties will be erased earlier which supports a gradual explusion of scalar 
gradients from small to large scales. In this vein, Fig. Q, much like the classic picture by 
Mclntyre, reproduced in Dritschel & Mclntyre [21], shows the inversion of a single "saw-tooth" 
using the general form in 0. As is evident — apart from the known asymmetry in the east- west 
jets in the beta plane equations — the more local the inversion (i.e. for smaller a) the narrower 
and stronger are the corresponding jets. 

To examine the nature of emergent flows for differing locality in the presence of dispersion, we 
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FIG. 4: Inversion of a single "saw-tooth" profile for varying a. As is evident, in addition to the east- west 
asymmetry, smaller a (more local) gives stronger and narrower jets. 

perform numerical simulations with spatially un-correlated initial data as shown in the first panel 
of Fig. Q. Setting e = 0.1, the scalar and zonal component of the velocity fields for a — 0.25, 1 
and a — 1.25 are shown in Fig. ([5]), quite clearly for a > 0.5 we have, what might be termed 
a coherent zonal flow. Note that, in accord with Fig. Q, the flows are broader and of smaller 
magnitude for increasing a. 

B. Resonant and near resonant interactions 

Evolving ([5]) in time, in the limit e —> 0, given the highly oscillatory nature of the integrals 
involved, the only terms that contribute to ^ are those which are exactly resonant. Following the 
notation of Chen et al. [25J, i.e. the first entry in a triplet (•, •, •) evolves via the interaction of the 
latter two entries, the limiting (e — > 0) dynamics are made up of : 

• Px — Qx — kx — i.e. the slow-slow-slow (5,5,5) interactions : All such interactions between 
purely zonal flows are clearly resonant. However, from Q we see that their interaction 
coefficients are zero. 

• k x = 0, p x , q x ^ i.e. slow-fast-fast (5, /, /) or fast-slow-fast (/, 5, /) interactions : It is easy 
to check that two modes with u ^ cannot feed (extract) energy to (from) one with uj — 0. 
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FIG. 5: 0,u fields for a = 0.25, 1 and 1.25 in the upper and lower panels respectively. In all cases e = 0.1. 
Note the finer scale of the flow as compared to the scalar field when a = 0.25. Further for increasing a, we 
obtain coherent zonal flows. 

In other words, two fast modes cannot generate a zonal flow. On the other hand, two fast 
modes can interact with each other using a zonal mode as a catalyst and this represents the 
passive advection of the fast modes via a zonal shear flow. 

• kxiPxiQx 7^ i.e. the fast-fast-fast (/,/,/) interactions : These are the three-fast wave 
interactions. Though they are algebraically difficult to achieve, these interactions represent 
fully 2D motion. 

As is evident, this situation is entirely analogous to the analysis of the beta plane equations as 
presented by Longuet-Higgins & Gill [26] . Further, there is a close correspondence between these 
resonant interactions and the limiting dynamics (i.e. when the Rossby number — > 0) of rotating 
three-dimensional (3D) flows. For the 3D rotating case [27], the analogous set of (5,5,5) inter- 
actions result in purely 2D flow (though now with a non-zero interaction coefficient); the second 
set of (/, 5, /) interactions leads to the passive advection of the vertical velocity via this 2D flow 
and inclusion of the third set of (/, /, /) resonances yields the so-called 2^ dimensional equations. 
Further, transfer from exclusively fast to slow modes in pure rotation is also precluded under 
resonance, i.e. (5, /, /) interaction has a zero interaction coefficient, where the zero frequency mode 
corresponds to a columnar structure. Therefore, in the limit e — ► (Rossby number — > 0) , it is not 
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possible generate a zonal flow (columnar structure) from exclusively non-zonal (non-columnar) 
initial data and forcing. 

Of course, the limit e —> is usually not of great physical interest, in fact we are usually more 
concerned with situations wherein < e < 1. This opens the door to non-resonant interactions, 
and quite naturally, we expect near-resonances — i.e. triads for which uj(p) + u(q) — u(k) = 5 <C e 
— to be of greater significance than the remainder that are, in terms of (5, e), far from resonance. 
As the (s, /, /) non-resonant triads have a non-zero interaction co-efficient we have the possibility 
of directly transferring energy from non-zonal to zonal modes. Choosing 5 for given e, such transfer 
to slow modes has been numerically verified in the beta plane and 3D rotating systems [28], [29J. 

IV. THE INVERSE AND FORWARD CASCADES 

With the aim of examining a systematic transfer of energy/enstrophy from one scale to another, 
we proceed to numerical simulations from weakly non-local to local conditions (i.e. a < 1) |12| , 
[30] . The numerical scheme is a de-aliased psuedospectral method with fourth order Runge-Kutta 
time stepping. All simulations are performed in Matlab at a resolution of 750 x 750. Further, we 
supplement the RHS of ([!]) with hyper-viscous dissipation of the form (— l) n+1 v n A n 9 (to act as 
a sink at small scales, with n = 4). As per Maltrud & Vallis [31J we choose v n = rms / (k m 2n ~ 2 ), 
where k m is the maximum resolved wavenumber. 

Denoting the domain, forcing and dissipation scales by k^^kf and kd respectively, we consider 
kL <C kf < kd (kz, < kf <C kd) which allows us achieve, for a given resolution, as large an inverse 
(forward) energy (enstrophy) transfer regime as is possible. The forcing is of random phase at every 
time-step and de-correlates with a time-scale r (this is chosen to be comparable to the fastest linear 
wave in the model). Further the forcing is localized to a few wavenumbers, and it is important to 
note — especially in the context of a small-scale case — that our forcing is not chosen in a "ring" 
of waveneumbers (i.e. it is not localized by means of exp{ — (\k\ — kf) 2 }) but rather, we strictly 
ensure that there is no projection of the forcing on small k x or k y (i.e. the localization is by means 
of exp{ — (k x — kf) 2 — (k y — kf) 2 }). This ensures that if we see the formation of zonal flows, it 
is due to the transfer of energy from modes that are, in a sense, "far from zonal" to zonal and 
near-zonal modes [32] . Finally, we also include a linear damping on the largest scales, of the form 
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— A#, which allows us to achieve a stationary state. In some of the inverse transfer simulations the 
time required to achieve stationarity is prohibitively large. In such cases the runs are halted when 
the power-law scaling ceases to evolve (i.e. only the very largest scales are still growing). 

A. Inverse transfer : Small-scale forcing 

As mentioned, in the context of the beta plane equations, one observes the spontaneous 
generation of zonal flows from isotropic random small-scale forcing (recent numerical simulations 
can be found in [22J, [22] )• In line with Rhines' argument for a non-isotropic streamfunction and 
dominant zonal flow, the reasoning — supported by numerical simulations [3T] . [35] — is that, 
the upscale transfer of energy is isotropic up to the so-called Rhines Scale, beyond which energy 
is preferentially deposited anisotropically in wavenumbers that satisfy k x /k y <C 1 [36j,[19j. In 
terms of the scaling, some forced simulations employing both large and small-scale sinks support 
the dimensionally predicted isotropic A: -5 / 3 and anisotropic A:" 5 power-laws for the non-zonal 
and zonal energy spectra respectively [37], |38| . Our aim is to examine the inverse transfer of 
energy for different a, and its sensitivity, for a given a, to the strength of the the dispersive term. 
Therefore we do not focus on the emergent zonal flows, but rather on the scaling associated with 
the inertial-range |39j under weakly- local to local energy transfer. 

Utilizing small-scale random forcing and large-scale damping as prescribed in the previous 
section, the energy spectra for different a with e = 0.5 are shown in Fig. Q. The non-dispersive 
estimate for an inverse cascade reads E{k) ~ fc - ( 7-2Q: )/ 3 [7J; quite clearly we obtain power-laws 
that are in reasonable agreement with the aforementioned dimensional estimate. Of course, our 
resolution is moderate and hence we do not pursue absolute quantitative accuracy in terms of 
ensemble simulations with appropriate error estimates. We now repeat these experiments with 
a dispersive term of differing strength (e = 1,0.1,0.05). The results for a — 0.5 are shown in 
the second panel of Fig. ([6]), as is evident, along with a shorter inertial-range the slopes become 
shallower with a progressively stronger dispersive term. The shortening of the inertial-range is 
expected as, it is easy to verify that the Rhines scale decreases with e. At present we do not have 
an explanation for the decrease in the slope of the power-law. Indeed, as the set of near-resonant 
interactions becomes smaller with decreasing e, it becomes prohibitively expensive to address these 
interactions in an adequate (or, at least, numerically consistent) manner. Hence, at our present 
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FIG. 6: The first panel shows energy spectra in the weakly-local to local inverse energy transfer regime 
with e = 0.5. We notice well-developed power-laws, and the slopes are in reasonable agreement with non- 
dispersive estimates. Note that these are slopes from individual realizations, not temporal or ensemble 
averages. The second panel shows the inverse transfer for a = 0.5, with e = 0.05,0.1,0.5 and 1. The 
sensitivity of the scaling with respect to e is evident, also note that the range up to which the scaling 
extends decreases with e. 

resolution of 750 x 750 we can only claim an agreement with non-dispersive estimates to be valid 
for e rsj 1. We also note the more pronounced steepening of the spectrum beyond the inertial-range 
for smaller e, though clearly we do not have adequate resolution to claim any spectral form for this 
steepening. 

B. Forward transfer : Large-scale forcing 

In the beta plane equations, the forward cascade of enstrophy has been studied by Maltrud & 
Vallis [3T] . It was noticed that the energy spectrum had non-universal properties, in particular 
they observed power-laws of the form k~ a with a E [3,3.5] [31] • From Q, it is evident that the 
dispersive term reduces the "effectiveness" of the nonlinearity; a conjectured consequence of which 
is the steepening of the energy spectrum from its dispersionless k~ 3 form [40J. Adopting a mixing 
perspective, the evolution of a large-scale initial condition (and its comparison with an identical 
passive field) in the beta plane equations was considered by Pierrehumbert [44]. Indeed, the rapid 
emergence of fine-scale features in the 9 field (implying a forward transfer of 9 variance) was seen 
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to yield a relatively high-wavenumber k~ x enstrophy spectrum. Much like our earlier discussion, 
the reasoning revolved around the passive driving of the high-wavenumber (or small-scale) features 
of the 9 field. Indeed, the k~ x spectrum coincides with the scaling of a small-scale passive field (in 
the convective range) driven by means of a large-scale smooth advecting flow, i.e. the Batchelor 
regime [H],[l6j. It is quite interesting that, the small-scale features generated from the initial 
large-scale 9 field do not spontaneously "re-combine" to, once again, result in the development of 
large-scale structures. In fact, the absence of any smooth, persistent coherent structures is also 
evident in the scaling, i.e. via the lack of anticipated spectral steepening that results from the 
presence such features. 

In addition to the increased physical locality seen in Fig. 0, from a spectral perspective the 
forward cascade is local for a < 1 [7], [47], (a = 1 is the well-known logarithmically divergent 
case [12], [IE]). In fact, the increased contribution from local triads (i.e. all legs of the triads are 
of comparable size) to the forward cascade has been noted in recent simulations comparing the 
SQG to the 2D Navier-Stokes case [49J. Possible consequences of the increased locality, such as 
the development of frontal discontinuities, the fractal nature of iso-9 level sets and the multifractal 
nature of the dissipation field have been examined in a series of SQG simulations [50] , [51] . As 
per non-dispersive estimates, the enstrophy spectra, shown in Fig. Q, have steeper slopes for 
increasing locality. Though, unlike the previous estimates for the inverse energy cascade, the 
slopes are clearly much steeper than the non-dispersive similarity hypothesis, which for a forward 
enstrophy cascade reads £(k) ~ fc _ ( 7_4a )/ 3 [7J. Repeating the experiments for e = 1,0.1,0.05 we 
always see the development of a clear inertial range; though now, as seen in the second panel of Fig. 
Q which shows the a = 0.5 case, the slopes are insensitive to e. Indeed, the principal difference is 
in the time required for a stationary state to emerge. Specifically, as the strength of the dispersive 
term increases (smaller e) the simulations have to be carried out for progressively longer durations 
(though, it is worth noting that these forward transfer runs are much faster than the previous set 
of simulations involving inverse energy transfer). 

V. CONCLUSIONS AND DISCUSSION 

We considered the evolution of a family of dispersive dynamically active scalar fields, the 
members of this family are characterized by the degree of locality in the — 9 relation. Further, 
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FIG. 7: Enstrophy spectra in the forward weakly- local to local enstrophy transfer regime for e = 0.5. We 
notice well-developed power-laws, but the slopes are significantly steeper than non-dispersive estimates. The 
second panel shows the same experiment for a = 0.5 but with varying e (e = 0.05, 0.1, 0.5 and 1 from top to 
bottom). Quite clearly, the slopes are fairly insensitive to e. Once again, we note that our estimates of the 
slopes are based on invidual realizations as opposed to a temporal or ensemble average that might be more 
suitable at higher resolutions for accurate quantitive slopes and error estimates. 

the dispersive term is strong, i.e. it is preceded by the inverse of a small parameter < e < 1. With 
regard to the effect of a on the energy /enstrophy transfer, according to the Fjortoft estimate, there 
is a monotonic relation between increasing a and the fraction of energy (enstrophy) transferred to 
large (small) scales. Further, the bias in energy /enstrophy partitioning increases with the spectral 
nonlocality of the transfer. As for the geometry of the scalar field, large isotropic coherent eddies 
are seen to develop from spatially un-correlated initial data for small a, while for large a we see 
the formation of a filamentary structure (much like a small-scale passive field when advected via a 
large-scale smooth flow). Rhines' argument for the spontaneous asymmetry in the streamfunction 
is seen to hold for all a > 0.5, while for a < 0.5 it is possible to maintain isotropy while satisfying 
the dual constraints of energy transfer to large scales and small frequencies. Indeed, simulations 
from un-correlated initial data confirm the emergence (and lack thereof) of dominant coherent 
zonal flows for a > 0.5 (a < 0.5) respectively [52J. 

Utilizing random small-scale forcing, for weakly non-local to local conditions (i.e. for a < 1), 
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we observe a clear power-law scaling in the energy spectrum, and for e ^ 1 the slope of the spectra 
agree with non-dispersive estimates. Under large-scale forcing, much like the inverse energy 
cascade, we see a clear forward cascade of enstrophy accompanied by a power-law scaling; though 
in contrast to the inverse regime, the slopes are much steeper than non-dispersive similarity 
hypotheses. Repeating the experiments with a dispersive term of differing strength shows that the 
scaling depends on e in the inverse transfer regime. Specifically, along with an expected shortening 
of the inertial-range, at our present resolution, the associated slopes also become progressively 
shallower as the dispersive terms becomes stronger. On the other hand, the forward transfer 
regime is quite insensitive to the strength of the dispersive term, in fact we observe the enstrophy 
cascade slopes to be fairly universal with respect to the range of e values considered. 

With regard to the mathematical aspects of ([I]), a = 0.5 is known to be special in the sense 
that it represents an open problem with regard to global regularity of non-dissipative solutions 
(see for example Constantin, Majda & Tabak [51] for an analogy between front formation in SQG 
and finite time singularities in the 3D Euler equations). In fact, present estimates for well behaved 
solutions require dissipation of the form A p with p — 0.5 for both, the non-dispersive and dispersive 
cases [55], [56]. Unfortunately, regularity results for 0.5 < a < 1 are presently un-settled. It is 
unclear if the difficulty in achieving regularity (by present techniques) arises abruptly at a = 0.5 
or whether one requires gradually stronger dissipation as a decreases from unity. Physically, this 
is interesting as a = 0.5 is precisely the border at which the advecting flow is of the same scale as 
the advected field. In fact, in the context of both the dispersive and non-dispersive cases, it would 
be quite interesting to know if a change in relative scales of the advecting and advected fields (as a 
crosses 0.5 from above) is connected to the deeper mathematical issue of the regularity of solutions. 
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